function j = getjac(a)
  j = [];
  